Repeated Bursts from Relativistic Scattering of Compact Objects in Galactic Nuclei 
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Galactic nuclei are densely populated by stellar mass compact objects such as black holes and 
neutron stars. Bound, highly eccentric binaries form as a result of gravitational wave (GW) losses 
during close fiybys between these objects. We study the evolution of these systems using 2.5 and 
3.5 order post-Newtonian equations of motion. The GW signal consists of many thousand repeated 
bursts (RB) for minutes to days (depending on the impact parameter and masses), followed by a 
powerful GW chirp and an eccentric merger. We show that a significant signal to noise ratio (SNR) 
accumulates already in the RB phase, corresponding to a detection limit around 200-300 Mpc and 
300-600 Mpc for Advanced LIGO for an average orientation BH/NS or BH/BH binary, respectively. 
The theoretical errors introduced by the inaccuracy of the PN templates are typically much less 
severe for the RB phase than in the following eccentric merger. The GW signal in the RB phase 
is broadband; we show that encounters involving intermediate mass black holes are detectable in 
multiple frequency bands coincidentally using LIGO and LISA. 



I. INTRODUCTION 

Close approaches between initially unbound compact 
objects (COs) can form bound binary systems if the 
gravitational wave (GW) emission, tidal dissipation, or 
interaction with other objects taps enough of the ini- 
tial kinetic energy [ll-Q. In particular, GW captures of 
black holes (BHs) and neutron stars (NSs) occur many 
times per Hubble time in dense stellar environments like 
galactic nuclei or globular clusters. These pairs provide 
sources for direct GW detection with Advaced LIGO or 
Virgo [H, [2 L and can lead to short-hard gamma ray bursts 
(SGRBs) Illil. 

The GW signals of these eccentric sources are very dif- 
ferent from standard quasi-circular inspirals. According 
to the leading order results 01 , for a given semimajor axis 
they are more luminous and are described by broadband 
spectra, which makes them detectable to larger distances 
and in a broader mass-range. Kocsis, Gaspar, & Marka 
(hereafter KGM) [l| investigated the detcctability of the 
GW burst emitted during a single passage, and found 
the SNR to be substantial only for encounters with a 
very small initial pericenter distance r^o ^ 6M (where 
M is the total mass in units G = c = 1), which oc- 
curs relatively rarely. O'Leary, Kocsis, & Loeb (hereafter 
OKL) included the much stronger GW signal from 
subsequent passages from bound systems following the 
first encounter, leading to an eccentric inspiral. Remark- 
ably, the expected detection rates of these sources for Ad- 
vanced LIGO is comparable to other types of waveforms, 
between 1-10'^ yr~^. The large uncertainty is mostly due 
to the unknown number and mass distribution of BHs in 
galactic nuclei (see Appendix [C] for further discussion). 

As the binary evolves from the initial very eccentric 
phase towards the less eccentric phase, the GW signal 
initially consists of well-separated repeated bursts (RBs) 
for minutes to days, and later transitions to a continu- 



ous inspiral waveform, a short but powerful chirp (OKL). 
The signal evolves from the RB to the chirp phase within 
the frequency band of Advanced-LIGO type instruments, 
making these GW signals particularly rich in features and 
very unique among other sources. These GW sources, if 
involving NSs, have electromagnetic counterparts, mak- 
ing them interesting candidates for multimcssenger as- 
tronomy [1, Q. 

Existing techniques are not well suited to dig these 
GW signals out of the noise, in either the RB or the final 
chirp phase. In the RB phase, individual GW bursts are 
relatively weak compared to the instrumental noise, mak- 
ing burst search algorithms insensitive to these sources. 
Nevertheless, since the time evolution of successive bursts 
can be predicted theoretically, and there are hundreds to 
thousands of well-separated bursts in the RB phase, it is 
in priniciple possible to optimize detection algorithms to 
coherently detect the full sequence of bursts. Regarding 
the final chirp, existing matched filtering searches with 
circular inspiral templates are also expected to be inef- 
fective, as here the eccentricity is still considerable [lol - 
|l3|. Post- Newtonian (PN) or effective one body (EOB) 
body waveforms have not been developed to sufficient ac- 
curacy for eccentric orbits of comparable- mass binaries 
with a small pericenter distance [l3l [l^ . Direct numer- 
ical experiments are restricted to non-extreme eccentric- 
ities (e < 0.7) and only a limited number of configura- 
tions have been tested [isl - fTsj . Without sufficiently ac- 
curate theoretical templates, matched filtering detection 
techniques will be prone to large theoretical errors p^ . 
These issues might be expected to be less severe in the 
RB phase, where the binary separation is relatively large, 
as long as the source is not in the zoom-whirl regime |20l - 

In this paper, we focus on the detectability of the GWs 
emitted in the RB phase. We examine the relativistic 
corrections to the evolution of the GW-capture binaries. 
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Wc numerically integrate the 2.5PN and 3.5PN equations 
of motion of the binary, including the radiation-reaction 
force. It is important to note that we do not include 
trajectories that technically whirl (execute a full 27r or 
more around periastron) since these orbits are by neces- 
sity in a regime where the PN equations of motion are 
unreliable. Numerical relativity is needed to examine ad- 
ditional boosts to SNR due to whirls. We calculate the 
GWs emitted and evaluate the numerical Fast Fourier 
Transform (FFT), to compare with the detection thresh- 
old of GW instruments. We study how the SNR accumu- 
lates in time during the evolution, and examine whether 
these broadband waveforms can be detected coinciden- 
tally in separate frequency channels with different GW 
instruments. (We find that they can.) We provide a 
brief estimate of the event rates in the Appendix. 
We use units G = c = 1 . 



II. EVOLUTION OF ORBITS 

We integrate the instananeous 2.5PN and 3.5PN equa- 
tions of motion of Will and collaborators, including spin 
corrections and dissipation due to gravitational radiation 
emission (26l - [30j . Simultaneously to numerical integra- 
tion of the trajectories, we calculate the two GW polar- 
izations. Note that this is different from the approach 
used for quasicircular orbits, where the orbit-averaged 
fluxes are calculated to a much higher order: 3PN order 
beyond the 2.5PN leading order flux [H [UHll] . Direct 
integration of the equations of motion allows greater flex- 
ibility when working with very high eccentricities, where 
averaging the GW flux over a Newtonian approximation 
to the orbital geometry would be very inaccurate. Q 

We use the 2.5PN and 3.5PN approximations to in- 
vestigate the relativistic corrections to the OKL study, 
and to assess the calculational uncertainties. While the 
3.5PN approximation is more accurate than the 2.5PN 
approximation at large separations, the PN calculation 
breaks down interior to rp < lOM in this calculation 
scheme, as in this case the 3.5PN correction dominates 
over the 2.5PN terms and leads to an artificial increase 
of the eccentricity [36| . The PN approximation is not 
well behaved in this important regime. However, the 
binding energy decreases monotonically until merger for 
the 2.5PN calculation. For the 3.5PN runs, we termi- 
nate the simulation where the magnitude of the 3.5PN 



^ Our approach is also different from that used in Damour et al. 
[sSl or Arun et al. |35| . which give the phasing of binaries to 
3.5PN order by averaging the radiation reaction over an orbital 
period and using this in a calculation of angular momentum flux. 
Orbit-averaged fluxes are computed to higher orders than ex- 
plicit equations of motion. As a result, variations in the signs 
of radiation-reaction terms at different orders can be washed out 
with orbital averaging. However, the limits of validity of the PN 
expansion are still pressed at close radial separations, regardless 
of the approach. 



perturbation terms first dominate over the lower-order 
terms. This is usually in the final chirp phase. Predic- 
tions for the RB phase are typically not affected by this 
truncation, as long as the binary is not in the zoom-whirl 
regime. 



A. Binary formation 

There are three possible outcomes after the first close 
encounter between two compact objects, depending on 
the initial conditions: unbound quasi-hyperbolic trajec- 
tory, capture into bound quasi-eccentric orbit, or direct 
collision. Here we examine relativistic corrections to the 
capture cross section. 

The event rates of these waveforms arc sensitive to 
the critical impact parameter for capture into a bound 
orbit 6niax = ^max('7^i, m2, w). Here mi and m2 are the 
component masses, and w is the initial relative velocity 
at infinity. To leading order, a bound (non-plunging) 
system forms if 
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where M = mi + m2. Here the upper bound assumes 
quadrupolar radiation emitted on a hyperbolic trajectory 
j37l |38| , and the lower bound is valid in the test- particle 
limit around a Schwarzschild BH (see KGM and OKL). 

Figure [1] shows the boundaries for 2.5PN and 3.5PN 
simulations for equal masses and for mass ratios 
toi/to2 = 0.1 and no spins. The green line shows that 
Eq. ^ used by OKL for binary formation is in excel- 
lent agreement with our numerical post-Newtonian cal- 
culation for these mass ratios. This is also expected 
from analytical orbital-averaged PN estimates for the 
typically nonrelativistic initial velocities in galactic nu- 
clei {w < 0.01)0 We are unable to resolve orbits that 
cross interior to lOM, even if they arc not direct cap- 
tures, as in this regime higher order PN effects may be 
more significant [36J. In the following wc focus on orbits 
that are captured in bound eccentric orbits. 

Combining the above estimates for the impact param- 
eter with the expected number density of objects in the 
galactic nucleus and their velocity distribution can be 



^ The post-Newtonian correction to the RHS of Eq. ([TJ was cal- 
culated by Junker and Schacfer [39j . Expanding their Eq. (48) 



in a power series in w to first to leading order gives 
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which yields deviations from the leading order term by less than 
3% for equal masses rj = 1/4 and w < 0.01. The correction is 
larger only for very unequal mass ratios, but such sources are not 
expected to exist for terrestrial GW instruments, based on the 
frequency limit of the instrument (implying M < 100 Mq) and 
the minimum mass of BHs and NSs (M > 1 Mq). 
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FIG. 1. Possible outcomes of the encounters depending on the initial velocity and impact parameter: escape (gray), capture 
into bound elliptic orbits (white), or direct collision (black). Lines show the analytical estimates based on Eq. ([1]). Top and 
bottom panels are for mass ratio 1 and 0.1, left and right panels correspond to 2.5PN and 3.5PN calculations, respectively. 
The 3.5PN calculations were terminated at radii where the derived perturbations are unphysical in the black region. 



used to make estimates on the likelihood of such encoun- 
ters. This excersize, summarized in Appendix [Cl yields 
that out of 10^ compact objects in a single galactic nu- 
cleus, only a few binaries form in a billion years. The 
corresponding instantaneous fraction of objects in the bi- 
nary forming region shown in Fig. [1] is extremely small 
on average. However, if the detectable distance of these 
sources is sufficiently large, the total rate from all observ- 
able galaxies may be quite high. 

We conclude that relativistic corrections do not mod- 
ify the capture cross section over the 10% uncertainty in 
the simple OKL estimate associated to the relative ve- 
locity distribution. These relativistic corrections are 
negligible compared to the theoretical uncertainties in 
the event rates as discussed in Appendix [Cj The detec- 



tion rates however may be affected by relativistic correc- 
tions through variations in the detectable distance of the 
source, which we investigate below. 



B. Orbital evolution 

Shortly after formation, the binary is very eccentric, 
the orbits are nearly radial. Due to GW losses near peri- 
center passage, the apastron decreases faster than the 
pericenter distance, leading to a decrease in eccentricity. 
The orbits exhibit large relativistic pericenter precession 
(see Levin, McWilliams & Contreras [s^ for a gallery of 
orbital trajectories). 

We compare our PN simulations with the leading or- 



0.001 0.01 0.1 . 1, 10 100 1000 

'.g[(V40M)^ (4n)-'M,„min] 




0.001 0.01 0.1 , 1, 10 100 1000 

f,g[(y40/M)' (4rir^M^„minl 



0.001 0.01 0.1 . 1, 10 100 1000 




0.001 0.01 0.1 . 1, 10 100 1000 



FIG. 2. The evolution of Vp and e as a function of time to merger, for mass ratio q = \, zero spin, using 2.5PN (left) and 3.5PN 
(right). The black lines are the analytic solutions of Peters [331 ■ Time is measured backwards from the innermost orbit of the 
simulation. Circles denote the transition from the RB to the final chirp phase for mi = 7712 = 10 Mq. 
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FIG. 3. The evolution of eccentricity as a function of pericenter distance, for mass ratio g = 1, zero spin, using 2.5PN (left) 
and 3.5PN (right). The black lines are the analytic solutions of Peters [37| . 
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der analytical formulas of Peters [l^l in Fig. [2] and |3l 
In OKL, it was demonstrated that for the later, assum- 
ing an initially parabolic orbit, e(rp/rpo) and eit/to) arc 
universal functions, independent of masses and the im- 
pact parameter. These physical parameters affect only 
the time (here t is measured from merger) and length 
scales (tt^^rpo). The evolution might be expected to be 
only slighly different for initially hyperbolic orbits, as 
long as the velocity at close approach is dominated by 
the gravitational binding energy, not the initial kinetic 
energy. The value of the relative velocity at infinity sets 
the maximum initial pericenter distance for binary for- 
mation. We examine the inaccuracies in the evolutionary 
curves due to relativistic corrections. 

Figure [2] and [3] plot the orbital evolution rp{t/to), 
e-it/to), and eirp/rpo) for 16 different initial {b,w) drawn 
from the white basin of Fig. [T] with w ~ 10~^c. The 
leading order analytic results arc shown as black curves, 
while the different panels correspond to 2.5PN and 3.5PN 
calculations. The figures show that both the 2.5PN and 
3.5PN calculation asymptote to the leading order eccen- 
tricity curve for large rp. For widely formed binaries, the 
deviations become significant interior to rp ^ 20M in the 
2.5PN calculation, or interior to rp ~ 40M for the 3.5PN 
calculation, while they are consistent to smaller separa- 
tions if the initial rpp is less. Part of the discrepancy 
between the simulated and the black analytic curves in 
Fig. [2] is that in the latter case time is measured from 
merger, while in the former it is measured from the last 
simulated orbit. This overall time shift is more prominent 
closer to merger on a logarithmic scale. Interestingly, the 
full 2.5PN calculation decreases the eccentricity steeper 
as a function of pericenter distance than in the leading 
order orbit-averaged approximation of OKL, while the 
3.5PN calculations is just the opposite, leading to a shal- 
lower eccentricity decrease. This is consistent with the 
results that radiation reaction is over-estimated at 2.5PN 
order and this over-estimate is tempered at 3.5PN or- 
der [36j . These two approximations bracket that used 
in OKL. For higher w, corresponding to the innermost 
regions of galactic nuclei, the highest initial rp values in 
the figures (dark blue curves) would not form binaries: 
Tpo^max ~ (93,50,25) for w = (0.001, 0.003, 0.01) c, re- 
spectively, (see Eq. (18) in OKL), but other curves with 
smaller r^o remain similar for different w. 

Circles in Figs. [2] and [3] mark the approximate bound- 
ary between the RB phase and the final chirp, where the 
time duration between individual GW bursts (i.e. the 
orbital time) At satisfies At > 5//min, where /min is 
the minimum frequency for a given detector. For Ad- 
vanced LIGO, /min ^ 10 Hz, so we require At > 0.5 s 
in the RB phase. Note, that from Kepler's law, this 
amounts to an approximate constraint on the semima- 
jor axis, a = rp/{l — e) > SIM^^^'^M, where M is the 
total binary mass, and Af20 = M/2OM0. 

Alternatively, we will also examine the signal to be 
in the RB phase if the GW signal is comprised of short 
duration bursts and longer silent periods, requiring that 



the silent periods are at least a given factor (e.g. 2 or 
4) larger than the burst duration. This later definition 
is equivalent to setting the eccentricity to be larger than 
some Cmin in the RB phase (e.g. emin = 0.45 or 0.6), 
independent of the semimajor axis. In this case the end 
of the RB phase is a horizontal line in Fig. |31 

Simulations with different mass ratios and spins lead to 
similar curves as those in Fig. [31 The 2.5PN and 3.5PN 
simulations with mass ratio g = 0.1 track the analyti- 
cal leading order curves for rp{e) to somewhat smaller 
pericenter distances (down to rp = (6, 7, 10, 15)M for 
fpo = (8,10,20,40), respectively). The orbit evolves 
to much smaller eccentricites before plunging. We have 
also run calculations with extremal spins in aligned, an- 
tialigned, and perpendicular configurations with respect 
to the orbital angular momentum. The result is qual- 
itatively very similar to the nonspinning case. In the 
aligned configuration, the time evolution of the eccentric- 
ity tracks the Newtonian result much more closely than 
in the nonspinning case, while the spin-orbit precession 
adds a slow periodic modulation to the evolution if the 
spins are initially perpendicularly oriented. 

The 3.5PN approximation leads to an apparent eccen- 
tricity increase at rp ~ lOM if e < 0.15. This feature 
is close to the point at which we truncate the simula- 
tion due to the breakdown of the approximation, so we 
take it with a grain of salt [s^ . Essentially we are seeing 
the apastron decrease more slowly than periastron rp 
leading to an increase in e = {ra — rp)/{ra -'r rp). For 
a dissipating orbit, eccentricity is not a precisely defined 
quantity, but the GW signal does show qualitative mea- 
sures of this e, such as a broadband character, noticeable 
for e > 0.1. In the extreme mass ratio case, eccentricity 
increase is known to occur only for marginally plunging 
orbits with rpo ~ 4 HI. The 2.5PN and 3.5PN calcu- 
lations agree until the end of the RB phase marked by 
circles in Figs. [5] and |31 The calculations are roughly 
consistent for the highly eccentric orbits e > 0.6, except 
for small r^o ^ 15 M which represent zoom-whirl orbits. 
In general, the true orbits are expected to whirl at small 
separations and exhibit larger precession than in the PN 
calculations. 



III. GRAVITATIONAL WAVES AND THEIR 
DETECTION 

We obtain the emitted GWs during the orbital evolu- 
tion as a function of time. In practice, we calculate the 
instantaneous /i+(t) and h^{t) polarizations of the strain 
amplitude in the direction of the orbital axis from the in- 
stantaneous phase space elements (for details, see [36|). 
For each mass ratio, we can use a single simulation to de- 
scribe sources with arbitrary total masses and source dis- 
tances, by scaling the amplitude and time proportionally 
with Mz/di^ and M^, respectively, where Mz = (1 + z)M 
is the cosmological redshifted total mass and c^l is the 
luminosity distance. 
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Given h^(t) and hx{t), the instrument measures a 
combination 



h{t) ^ F+h+{t) + F^h^it) 



(3) 



where Fj^. and Fx are the antenna beam pattern coeffi- 
cients, which depend on the orientation of the detector 
with respect to the binary (see Eq. (104) in Ref. and 
satisfy < |-F+,x I < 1, where F+ = 1 if Fx =0 (and vice 
versa), and on average {F^ x) = 1/5- We discuss how we 
make inferences for an average binary orientation from 
a waveform corresponding to the optimal orientation in 
Appendix O 

The simulations confirm the qualitative waveform fea- 
tures presented in OKL. During the first part of the evo- 
lution, when the eccentricity is very large, the binary 
emits repeated GW bursts (RB) during successive close 
approaches. The relative amplitude of the two polariza- 
tions are modulated by GR precession. The time sepa- 
ration between successive bursts decreases rapidly on a 
logarithmic scale, as the eccentricity decreases. The time 
duration of individual bursts changes much more slowly. 
Eventually, as the eccentricity becomes small, the RB 
phase evolves toward a continuous chirp signal. 

In order to assess the detectability of the waveforms, 
we calculate the numerical EFT of the waveform for both 
polarizations, /i_|_.x(/), and compare to the sensitivity 
level of GW instruments. This is more accurate than 
the OKL estimate, as that relied on the stationary phase 
approximation to estimate the Eourier amplitude of each 
harmonic. 

The SNR is given by 
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where Sn{f) is the one sided spectral noise density and 
h{f) is the sky position, binary orientation, and polar- 
ization averaged GW signal spectral amplitude (see Ap- 
pendix 

Eq. (|4]) shows that 2fh and \/ fSn correspond to the 
angular averaged spectral signal amplitude and RMS 
noise amplitude per logarithmic frequency bin, whose 
ratio gives the SNR per logarithmic frequency bin. In 
the following we compare these dimensionless quantities 
when discussing the detectability of the signal. 
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FIG. 4. The characteristic spectral amplitude for circular or- 
bit (zero spin) initial conditions ro = 80Af, for 2.5PN (blue) 
and 3.5PN (green) and masses mi = m2 = 10 M0 (top 
curves) and mi = 0.1m2 = 1.4 M© (bottom curves). Dot- 
ted lines show the analytical stationary phase approximation 
using the 2.5PN orbit-averaged flux of Poisson & Will [ZJ. 



asymptote the orbit-averaged flux spectra for large sep- 
arations or small frequencies, but at higher frequencies, 
they lead to systematically lower and higher GW spec- 
tral amplitudes, respectively. Interestingly, the orbit- 
averaged flux result is well off our 2.5PN spectral am- 
plitude at / < 50 Hz (i.e. r < 20M for BH/NS bi- 
naries), but it is very close to our 3.5PN calculation. 
The very strong final peak and the final upper har- 
monics are artifacts of the 3.5PN calculation, as the 
orbital evolution slows down there around lOM, where 
/cut = 10-3/2;r-iM-i = 100 Hz X {M/20Mq)-^ and 
the eccentricity increases in the calculation (see Fig. [3|) . 
A similar spectral increase was identified by Buonanno 
et al. [13 using a different approach, where the GW fiux 
was calculated to IPN beyond leading order, correspond- 
ing to 3.5PN in our calculation, but it is not present using 
higher order corrections to the flux. Such higher order 
orbital averaged flux calculations are consistent with nu- 
merical simulations for quasicircular inspirals [H, l3ll - [33j . 
In the following we present both the 2.5PN and 3.5PN 
calculation results for eccentric orbits to gauge the error 
of our calculations in the RB phase. 



A. Circular orbits 

We first examine the convergence of 2.5PN and 3.5PN 
calculations for circular initial conditions. Figure S] shows 
the GW spectra in this case, upper and lower curves 
correspond to BH/BH and BH/NS binaries. The dot- 
ted lines show the analytical spectra for 2.5PN orbital- 
averaged flux calculation for circular orbits with the sta- 
tionary phase approximation (SPA), see Eq. (jA3[) in Ap- 
pendixIXlbelow. Both the 2.5PN and 3.5PN calculations 



B. Time evolution of the GW spectra 

Figures [5] and M show the angle-averaged GW spec- 
tral amplitude, for different initial pericenter distances 
fpo = (17, 33, 55)M, using (2.5, 3.5)PN orders, and mass 
ratios q = (1,0.1). The colored curves represent sep- 
arate minute and hour segments as marked, while the 
top green curve shows the total root-sum-squared (RSS) 
spectra of the full waveform. The dotted line shows the 
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FIG. 5. The characteristic spectral amplitude for GW capture orbits for 2.5PN (panels in row 1) and 3.5PN (row 2) calculations, 
and masses mi = m2 = 10 Mq for rpo = (17, 33,55)M (left, middle, and right panels) for a binary at 100 Mpc with an average 
orientation and zero spin. Top green curves correspond to the full waveform, others show the contribution of 1 minute (left 
and middle) and 1 hour segments (right panels), respectively. Thick red line is the Advanced LIGO sensitivity. 




FIG. 6. Same as Fig.[5]but for mi — O.hno = 1.4M0 resembling BH/NS encounters. 
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spectral amplitude for circular inspirals for reference (see 
Eq. (|X3)) ). The signal is in the Advanced LIGO band 
and well above the sensitivity level for minutes to hours. 
Initially, in the RB phase, the waveforms are broadband 
within a short-duration followed by long silent periods. 
Although the orbital frequency /orb is well outside the 
LIGO band hcrcH, the characteristic frequency of the 
bursts are the inverse timescale of pcricenter passage, 
fp - /orb(l + e)i/V(l - e)3/2 » /orb. Remarkably, the 
signal transitions from the RB phase to a chirp signal in 
the detector band, as the eccentricity quickly decreases 
in time. 

Next, we split the GW signals into approximately 1 
second segments (see Appendix |B] for details) and cal- 
culate the angular-averaged SNR for Advanced LIGO 
in each segment. Figure [7] shows how the SNR accumu- 
lates when measuring the signal from the first passage to 
a time before merger for binaries at 100 Mpc. Dif- 
ferent curves correspond to binaries with different peri- 
center distances at first passage rpQ as labelled (same as 
in Figs. [5] and Different panels correspond to dif- 
ferent binary masses and FN order (see figure caption). 
Initially, in the RB phase of the binary evolution, the 
SNR accumulates in discrete bursts. Although hard to 
see on the logarithmic scale, the strength of successive 
bursts is nearly equal. Individual bursts are detectable 
separately at high significance only if the first passage is 
sufficiently close (rpo < 15 for a sky-averaged S/N > 5 
for mi = 7712 = 10 Mq at 100 Mpc). The trends and 
order of magnitudes arc broadly consistent with KGM. 
The end point of the RB phase is marked with big cir- 
cles, where the GW signal becomes continuous, starting 
the final chirp. 

The figure shows that the total SNR of the RB phase 
can be a substantial fraction of the total SNR. This pre- 
diction is robust in both 2.5FN and 3.5PN calculations. 
Therefore we conclude that a coherent search for the train 
of bursts has a potential for detection even if the SNR of 
individual bursts is small^ 

C. Total signal to noise ratio 

Up to this point we restricted to a single choice of 
masses for BH/BH and BH/NS mergers. Now we extend 
the analysis to all possible total masses observable to 
LIGO between 2M0 to 1000 Mq. 

The contours in Figures [5] and O show the angular- 
averaged SNR for Advanced LIGO in the RB phase for 



^ In fact, depending on mass, fp and /oj-b can be in the LIGO and 
LISA band coincidentally in the RB phase, see Sec. IIII El below. 

* Note that this implies that the RB phase does not show up clearly 
on time-frequency plots of the SNR, since the signal power is 
small during each passage. Such time-frequency plots are more 
useful for individual burst searches or for the final chirp, when 
the total signal duration is not very large, and the SNR per unit 
time is substantial. 



sources at 100 Mpc, as a function of r^o and M, for mass 
ratio q = 1 and 0.1 (top and bottom panels), for 2.5PN 
and 3.5PN calculations (left and right panels). In Fig- 
ures [5] and [21 we select the segment of the signal where 
e > 0.6 and the e > 0.45, respectively, and evaluate the 
corresponding numerical FFT and SNR for Advanced 
LIGO. Note that the pericenter passage timescales in 
these cases are at least 5 times and 3 times shorter than 
the orbital timescale, so the GW signal consists of well 
separated RBs between silent periodsQ The SNR re- 
sults for the 2.5PN and 3.5PN calculations in the RB 
phase are consistent to within 30%. The 3.5PN curves 
cannot resolve the SNR when r^o ^ 10 for reasons al- 
ready mentioned (namely, the poorly behaved approxi- 
mation). We expect the SNR to be significantly larger 
for zoom- whirl orbits with Tpo ^ lOM relative to our cal- 
culations (KGL). As shown by Fig.lH the 2.5PN calcula- 
tion systematically undcrpredicts the SNR at all frequen- 
cies by a factor ^ 1-2 for circular orbits at frequencies 
/ ^ 0.3/iscOj while it may be even more uncertain at 
higher frequenciesc closer to ISCO. 

Comparing Fig.[8]to Fig. 11 of OKL, we see that a con- 
siderable fraction of the total SNR is in the RB phase, 
for masses less than ^ 20 Mq -|- 20 Mq and initial peri- 
center distances rpo < 40, making the RB signals typ- 
ically detectable to several 100 Mpc with high signifi- 
cance. Our highest SNR results in the RB phase, are 
typically a factor 2-3 lower than OKL's for the full wave- 
forms. While OKL claimed that intermediate mass BHs 
with mi = 1712 = 400 Mq are detectable to 1 Gpc with 
SNR = 5 with Advanced LIGO for r^o ~ 6M, our 2.5PN 
and 3.5PN calculations cannot accurately model orbits 
in this range to either confirm or rule out such claims. 



D. Precession effects 

Figures [SHZ] assume orientation and polarization aver- 
aged waveforms for zero spin. In reality, however, each 
GW detector will be sensitive to a single linear combina- 
tion of the -I- and x polarizations. The measured GW 
signals are strongly modulated by the IPN GR precession 
of the eccentric orbit within the orbital plane. Further, 
if the objects are spinning with general non-aligned spin 
orientations, spin-orbit precession further modulates the 
signal waveform at 1.5PN order [i^. In the leading order 
flux averaged approximations, the precession periods are 



Note however, that the time duration between individual bursts 
may be smaller than 0.5 sec in this case, so the circles in [2l [3l 
and[7]do not coincide with e = 0.6 or e = 0.45, see Fig. [3] 
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FIG. 7. Angular-averaged cumulative SNR for Advanced LIGO from the first passage as a function of time to merger for 
2.5PN (left) and 3.5PN (right) calculations and binary masses mi = m2 = 10 (top) and mi = 0.1m2 — 1.4 M© (bottom). 
Different curves correspond to different initial pericenter distances as labelled, same as in Figs. [2] and [S] Far from the merger, 
the binary is in the RB phase, and the SNR accumulates during close approaches. Circles resemble the transition from the 
RB to a continuous chirp signal, where the orbital time is 0.5 s. All panels assume a source at 100 Mpc, and average binary 
orientation. 



respectively, In contrast, the total signal duration after the first fly by 

is approximately (see Eqs. (13) and (27) in OKL) 



2t: r', 



"1^ ^ - " = 0.89 hr X (4r;)-3/2 (Jj^S"'" . (g) 



MV-^ (1 - e)3/2 

M / r„ \5/2 1 + 



20 Mq V30M/ 

^'^^ ^ 2Q y[ i 30 M ) 7l — g\3/2 ' The Earth spin also modulates the measured waveform, 

which could be significant for waveforms lasting sev- 
^ _^r^ n + e V'^ eral hours, i.e. if rpo > 30M X (A//20Mo)-4/2i (477)2/7. 

yi _ Clearly, the evolution time in the RB phase is almost al- 

M 3/1 \ ways much larger than the precession timescales and 

= 8_4s X ( ^ I ) ■ (5) ^^^^ Fig. [71). Accounting for these effects is going to 

20 M0 \30M/ \l — ej be crucial for detection algorithms. 
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FIG. 8. SNR for Advanced LIGO in the RB phase when e > 0.6, as a function of the pericenter distance of the first approach 
rpo and binary mass for mass ratio q — I (left) and 0.1 (right panels), in the 2.5PN (top panels) and 3.5PN calculation (bottom 
panels). The source distance is 100 Mpc. For other distances the SNR is reduced proportionally. 



E. Coincident multiwavelength observations 

An interesting unique characteristic of these GW sig- 
nals is that they arc broadband, and can be detected 
with multiple GW detectors coincidentally in different 
frequency bands. Figure [TU] shows the GW spectra for an 
optimally oriented (i.e. face-on) binary for total masses 
200 Mq < M < 720 Mq (intermediate mass black holes 
IMBH) and 7x lO^M© < M < 3x IO^Mq (supermassive 
black holes SMBH) with different impact parameters at 
100 Mpc. For SMBHs, the maximum observation time is 
limited to the final 10 years before coalescence. The fig- 
ure shows that the spectral range of the signal spans the 
frequency range of multiple instruments for these masses 
for a wide range of initial Vpo- Ultimately, the frequency 
range of the Einstein Telescope is most ideal to detect 
these coincident LIGO/LISA sources with IMBHs. Go- 
ing to even higher masses and lower frequencies, we find 
that parabolic GW captures of SMBHs are detectable 
with future pulsar timing arrays (PTAs) such as the SKA 



and LISA, coincidentally. 

The event rates of such IMBH or SMBH encoun- 
ters within a few tens of Mpc is currently unknown. 
Portegics-Zwart ct al. |4J] predict that the galactic cen- 



ters hosts 50 IMBHs, which could undergo GW capture 
events. GW observations of these events could prove the 
existence of such a population. Regarding SMBHs, they 
are much more luminous and easier to identify using elec- 
tromagnetic (EM) observation. There is an ongoing effort 
to search for SMBH binaries using EM observations. 3C 
66B, a nearby radio galaxy at 80 Mpc, was interpreted 
as an SMBH binary with a total mass 5.4 x 10^° M© 
[ist . However, Jenet et al. [1^ have shown that the 
GWs would be detectable with PTAs but are not ob- 
served, ruling out the SMBH binary interpretation for 
this particular source. Based on our SNR estimates, we 
conclude that PTAs could search for RB sources at sim- 
ilar distances. 

The prospects for multiwavelength observations may 
be somewhat better in reality than shown in Fig. 1101 Our 
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FIG. 9. Same as Fig. [8]but for e > 0.45, corresponding to fp/fovh < 1/3. 



calculations are artifically truncated at relatively small 
frequencies / < 0.4/isco due to the inaccuracy of the 
simulations at small seperations r < lOM. The com- 
plete signal, including the final chirp and the following 
ringdown, extends to higher frequencies than shown in 
Fig. [TUl and extending into the LIGO and LISA bands 
for IMBH and SMBH binaries, respectively. As the in- 
spiral and ringdown phases are both detectable, these 
sources constitute a new class of "golden binaries" , and 
may be useful to probe strong field gravity (47j 



IV. DISCUSSION 

We have examined the evolution and GW spectra of ec- 
centric binaries formed by GW emission using the 2.5PN 
and 3.5PN equations of motion of Will and collabora- 
tors [26l - [30l ISq . The capture cross section, evolutionary 
tracks, and the GW spectra are in remarkable agreement 
with the simple analytic estimates of OKL. The 2.5PN 
and 3.5PN results bracket those in OKL. 

After the formation of the binary, the GW signal is 



described by a long repeated burst (RB) phase lasting 
minutes to days, followed by a continuous powerful chirp. 
The signal evolves from the RB to the chirp phase within 
the frequency band of Advanced-LIGO type instruments. 
The maximum distance of detection for waveforms in the 
RB phase is around 300-600 Mpc for Advanced LIGO for 
an average orientation BH/BH binary with SNR ~ 5-10 
for 10 M < Vpo < 25 M and 20 Mg < M < 100 Mq (see 
Figs. [5] and HI). We find that the total SNR is substantial 
already in the RB phase when the eccentricity and sepa- 
ration are relatively large. Numerical relativity may best 
resolve the prospects for detecting the GW signal in the 
final powerful chirp phase, or signals with smaller impact 
parameters leading to zoom-whirl orbits. 

We found that rclativistic corrections do not greatly 
modify the event rate estimates of these waveforms. OKL 
have shown that the Advanced LIGO detection rates for 
these sources may be around 1-3000 yr^^, depending on 
the number of stellar mass BHs with masses between 10- 
40 Mg in galactic nuclei. The event rates for IMBH en- 
counters and encounters involving a BH and NS may be 
equally numerous. Discarding the final powerful chirp. 
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FIG. 10. The GW spectra for supermassive and intermediate- 
mass BH binaries. The signal is marginaUy in range for coin- 
cident multiwavelength GW detections with LIGO, LISA, and 
SKA. Different overlapping curves represent different impact 
parameters corresponding to Figs. [3] and [S] The maximum 
observation time is limited to 10 years. 



and assuming conservatively that sources are detectable 
in the RB phase to 1 /3rd of the distance of the full sig- 
nal, the event rates for these sources may be reduced by 
a factor 27, to give a total detection rate 7?,rb ~ 0.03 
- 100 yr^^. The event rates in reality may be much 
higher than these estimates which in OKL was based on 
an isotropic density distribution of compact objects in 
galactic nuclei. However, heavier objects segregate into 
very anisotropic configurations efficiently through vector 
resonant relaxation |48| . Indeed, the observed distribu- 
tion of massive stars in the Galactic center is anisotropic, 
comprising two disks with a thickness of 10° [i^. This 
thicness is consistent with the prediction of mass segre- 
gation for objects with these masses in statistical equilib- 
rium with lighter 1 Mq stars. It is plausible to expect a 
similar anisotropic distribution of BHs in galactic nuclei. 
Since the event rates arc proportional to the squared den- 
sity of compact objects, one might expect that the true 
event rates are larger by up to a factor 100 than in OKL, 
making these sources much more numerous than other 
LIGO- VIRGO sources js^l (Kocsis & Tremaine in prepa- 
ration) . 

We have shown that encounters between IMBHs may 
be detected with Advanced LIGO and LISA coincidcn- 
tally if the source is within 50 Mpc. Similary for SMBHs, 
coincident detections may be possible with LISA and fu- 
ture Pulsar Timing Arrays such as the SKA. 

We have also estimated the SNR for detecting the GWs 
from BH/NS encounters, and found that a detection with 
SNR = 10 may be possible to 300 Mpc. These encounters 
may lead to tidal disruption events and may exhibit lu- 
minous coincident electromagnetic variations Q. Indeed, 



short-hard gamma ray bursts (GRBs) arc modelled as the 
merger of two NSs or a NS into a stellar mass BH, fol- 
lowing a circular inspiral. It is plausible to expect that 
GW captures leading to eccentric coalescences also result 
in similar phenomena. An important difference for very 
eccentric orbits is that the pericentcr separation may be 
smaller than the ISCO for circular orbits [5|. This may 
lead to tidal stripping, partial disruption, or the shatter- 
ing of the NS crust during close approaches [1, Q . As the 
mass of the NS is reduced, its radial size increases, so that 
tidal stripping becomes more efficient for successive close 
approaches. GRBs with observed precursors [5l| might 
correspond to these eccentric events. Thus, this process 
has a potential to generate electromagnetic bursts track- 
ing the GW signal. The LIGO data near GRBs could be 
searched for these particular GW signals. If such coun- 
terparts are successfully identified, such processes could 
be used as standard sirens to constrain the cosmological 
model [52, 53], and the mass of the gravitonin alternative 
theories of gravity [13] . 

Our 2.5PN and 3.5PN calculations demonstrate the 
slow convergence of the PN expansion for these encoun- 
ters in the final inspiral phase. More accurate calcula- 
tions would be necessary to make more accurate predic- 
tions on the detectabiltiy of these signals. Indeed, the 
waveform modelling precision may be an important lim- 
iting factor for concrete detection techniques and param- 
eter measurement accuracy. The theoretical errors due to 
the imperfect modelling of these signals may be dramatic 
for sources with small initial pericenter distances or dur- 
ing the later parts of the signals approaching merger [l^ . 
However, we have shown that the initial RB phase of the 
GW signal carries a considerable total SNR. While the 
theoretical modeling of the GW signal in this phase may 
be more accurate, their detections requires searching for 
a train of GW bursts over long timescales which have in- 
dividually a small amplitude. The standard LIGO- Virgo 
detection pipeline is not sensitive to these signals. A 
long-duration transient search with a network of instru- 
ments might be a more promising avenue for detection 
[ssj . However, as the expected waveforms are well de- 
scribed in the RB phase, optimized data analysis tech- 
niques could be developed for their targeted detection. 

Future studies should investigate how accurately can a 
network of detectors measure the physical parameters of 
these sources. These are the component masses, binary 
distance from the Earth, binary orientation, sky loca- 
tion, time and phase at merger, initial pericenter distance 
(or cquivalently, the impact parameter), and the initial 
velocity before the first encounter. The later two pa- 
rameters {rpo,w) affect the eccentricity evolution, which 
is nonnegligible when the signal is in the detector's fre- 
quency band. Measuring the eccentricity evolution yields 
Tpo, and is mostly insensitive to w. The later might be 
hard to detect directly, unless the GW signals from the 
first few passages can be resolved. However, since binary 
capture implies a maximum r^o for a fixed w, an estimate 
of fpo puts an upper bound on w. This may already be 
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sufficient to distinguish between sources in galactic nu- 
clei where typically v > 1000 km/s from those in globular 
clusters where v < 60 km/s (OKL). Regarding sky local- 
ization, the modulation related to GR and spin precession 
or the rotation of Earth, can be a substantial help, as the 
GW signal amplitude changes significantly due to these 
effects during these long duration signals. We expect this 
to greatly improve the measurement accuracy for these 
eccentric sources beyond the 10 deg^ accuracy of regular 
circular inspirals [56[. 
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Appendix A: Angular average waveforms 

We calculate the numerical waveforms in a direction 
along the orbital axis of the binary, and use it for a 
proxy for general inclinations. In particular, in the cir- 
cular orbit-averaged case, the radiation is dominated by 
the {i,rn) ~ (2,2) spin- 2 weighted spherical harmonic, 
so that h+{t) = (1 + cos^ t)Asin(27r/i) and hx{t) = 
(2cost)Acos(27r/t) as a function of the inclination angle 
relative to the line of sight. Here A gives the amplitude 
scale, which in the optimal orientation cos l. — 1 satisfies 
SA^ = {h"^ + /i^)opt- Averaging over an isotropic dis- 
tribution, 2{hl) = (28/15)^2 and 2{hl) = (4/3)A^ so 
that (hl+hl) = {2 / 5) {hi + hi) opt; the GW power in the 
optimal orientation is 2.5 times larger than the average 
GW power. We find that these identities are satisfied 
to within 1% accuracy in the Peters-Matthews leading 
order orbit averaged flux calculation for arbitrary cxccn- 
tricity. We will use them to change from the calculation 
evaluated at the optimal orientation to estimate the GW 
signal in the "average" case. 

The angular-average waveform amplitude is then 



h^if) = (Flhlif) + Flhlif)) = ^[h\U) + hl{f)Ut 

(Al) 

where we have used {F"^) = (^x) = 1/5 for the beam 
patterns and a factor 2/5 for the inclination-averaged 
GW amplitude. While in general, {Flh\) ^ {Fl){h\), 



however, the error introduced by this approximation is 
only 6% in the circular case, so we shall use Eq. (|Aip as 
an approximation. 

Note, that the signal measured by a single detector in 
the optimal orientation with F^ = 1 and F-^ = is larger 
than h{f ) by 

In the stationary phase approximation, the angular av- 
eraged spectral amplitude is 

'« = ^(l) '-'=^r"' (A3) 

assuming circular orbits, leading PN order orbit- 
averaged flux, and no spins [4l|. Here, the 4/5 fac- 
tor is the RMS of the polarization amplitude Ap = 
(1 -I- cos^ t)F_|_ — 2icos6Fx for random orientations. A 
similar approximate formula is available for eccentric or- 
bits summing over the contributions of orbital harmonics 
(OKL). 



Appendix B: Subtleties with numerical FFT 

There are some numerical subtleties we have to con- 
sider when calculating the FFT of a very long data 
stream, of length T, spanning several months in some 
cases with several kHz sampling frequency. To save com- 
putation time, we split the waveform into segments of 
size AT. side effects: (i) the spectrum becomes inacces- 
sible at / < 2/ AT, (ii) the frequency resolution will be 
discrete with a fundamental frequency 1/AT, and (iii) 
it can introduce artificial features and numerical errors if 
the signal and its derivatives do not vanish at the bound- 
aries. To minimize these errors but optimize the compu- 
tation time, we apply the following procedure: 

1. Keep the segment size as large as possible. In 
practice AT, - AT = 2 x lO'^sec(M/2OM0) typi- 
cally for the i"^ segment, except for Fig. [3 where 
AT, - AT = 1 sec. 

2. Choose split points measuring them from the end 
of the timeseries (i.e. near merger), so that the last 
split interval is at least size AT. 

3. Adjust the split points between intervals to the 
nearest local minima of /i+(t) -I- h1,{t). 

4. Apply a gradual fade in and fade out near the edge 
of the waveform over a timescale ifadc ^ 500M. 

5. Resample the /i+(i) and h^ it) time series with uni- 
form time steps. This is necessary since the simu- 
lations use adaptive time steps. 

6. If the orbital time torb is larger than AT, then trun- 
cate the data stream to the periccnter passage, cen- 
tering the split window of duration T there. 
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7. Append the data stream with = hx{t) = 0, 
if the total signal duration T is smaller than AT. 

8. Calculate the FFT of the resulting time scries, cor- 
responding to h^{t) and hx{t), for each time seg- 
ment. 



Appendix C: Event rates 

Here we provide a simple estimate of the event rates of 
these sources in galactic nuclei, and highlight the main 
sources of uncertainty. For a more detailed treatment see 
O'Leary et al. 0. 

The mean number density of objects with mass m, in 
a dynamically relaxed galactic nucleus around a super- 
massive black hole (SMBH) of mass M, is 



3 - a, 2k,M, _3 

Mr)) = — ^ 

Att to, \ r,- 



(CI) 



where the total mass of stars within the radius of influ- 
ence, ri, is 2M,, the total mass fraction in to, is k,, and 



GM, GMo f M, 

Mo 



l-(2/fc) 



(C2) 



where we have used M, = Mo((t/(To)*^ according to the 
M, — a relation, in which ctq = 200km/s, Mq = 1.3 x 
10^ Mq, and fc = 4 [13, El- We extrapolate this density 
profile inwards until a radius which encloses only 1 BH, 

r.„,in = r,i = A^,-^/('-"-V,, (C3) 



where N, = 2k,A'/,/to, is the number of BHs within r.j. 

The maximum impact parameter, ^max for a GW cap- 
ture for given relative velocity w is from Eq. (17) of 
O'Leary et al. 



(C4) 



where rj is the symmetric mass ratio and TO.tot is the total 
binary mass. We will take equal mass binaries for which 
77 = 1/4 and m,tot = 2to, . In the following we assume for 
simplicity that the relative velocity is the circular velocity 
at radius r from the SMBH, w = v{r) = (M./r)^/^^ 

The scattering cross section is a^s = '""^max- The binary 
capture rate in a spherical shell of thickness dr/r is 



din?' 



F = (47rr^)<crcsW 



FoK, 



Ml. 

Mn 



(53/14)-2q. 



, (C5) 



in the second equation we have plugged in Eqs. (jC2HC4p . 
and 



\ c 

6.13 X IQ-^yv-^ X (3 - a.)^^ 



where ^ = n1{ri) / {n,{ri))'^ . FromEq. (|C5[) it is clear that 
for fixed r/r^, the rates are independent of m,, if their to- 
tal mass fraction in the cluster, k,, is fixed. This follows 
from the fact that dF/dlnr oc 7^^&,p^,x where (n,)'^ cx jti^^ 



and 6^ax see Eqs. (|C1|) and (|C4p . 

Integrate dF/dr between r, min ^ r.i to get the total 
rate in one galaxy, using Eq. (|C3[) . 
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where F^, = [2a, - 53/14] -^Fq and 
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(C8) 

Let us introduce normalized parameters with numbers 
representative of a Milky- Way-size galaxy, R, = k/0.025, 

M, = M4e6 X 4 X IO^Mq, to, = TO,i X 10 M0, r,,nin = 

n/N,. 

FiGN « 3.0 X 10-9 yr-i 60 ^1'/'' to."//" (C9) 

where ^ = 30 ^30, and we assumed a, = 2. The rate per 
Milky-Way galaxy (k, = M4e6 = to,i = 1) is around 
3 X 10-9^0 yr-^ 

The total rate is the sum of the rates of individual 
galaxies within a detectable distance, dmax, 



47r 



•^fnax^airiGN 



. -1 ^ -31/14 - 

= 4.5 yr ^30 k, to. 



M^/f ngal,5<ax,2 (CIO) 



where (iniax2 = rfmax,2 X 2 GpC, and Tlgal = ngal,5 X 

0.05 Mpc-^ This is comparable to the rates found in 
Table 1 of O'Leary et al. 

The take-away from these calculations can be summa- 
rized as follows. 

• In O'Leary et al. 0, Fig. 11, shows that the max- 
imum detection limit of the full signal (including 
the RB phase and the final chirp) is between 1- 
3Gpc for 10 M0 < to, < 500 Mq for a broad 
range of impact parameters for Advanced LIGO 
with S/N = 5. Eq. (|C9p shows that the rates are 
weakly sensitive to to, as long as their total mass 
in the cluster is fixed. The event rates involving 
a few IMBHs may be equally numereous as rates 
among many stellar mass BHs with the same total 
mass. 



^ To directly compare the rates in Table 1 of O'Leary et al. 
the values should be scaled up by a factor of ~ 2 to account for 
the slightly smaller normalization used in that paper. 
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• O'Leary et al. [2] have shown that the event rates 
are dominated by close first encounters, where the 
signal enters the LIGO band in the RB phase. This 
is to be expected as the event rates are dominated 
by the innermost objects, where the velocity disper- 
sion is large, requiring a close approach for binary 
formation. 

• Eq. (|C9[) shows that the rates are weakly sensitive 
to M,. This is in stark contrast to other gravita- 
tional wave sources whose rates scale proportion- 
ally to the mass of the galaxy. The net rates are 
determined by the number density of lower mass 
galaxies which greatly outnumber Milky- Way sized 
galaxies. 

• Observations show that many galaxies can signifi- 
cantly deviate from the M — a relation, in terms 
of n(ri). This implies that the average value of 
(0 = {nl)/{n.f > 1 in Eqs. ^ and JCH). If 
the RMS scatter in the densities from the M — a 
relation is 5, then ^ = 25. 

• The exact average value of and a, is uncer- 
tain, due to uncertainties in the final initial mass 
function in these environments. While k, = 2.5% 
(k = 1) may be reasonable for the center of the 
Milky Way, assuming 20, 000 BHs of mass 10 Mq 
(see refereces in 0]), other values may also be pos- 
sible in general. Indeed, the relaxation timescale 



in M32, which is a dwarf elliptical galaxy hosting 
a 2 X 10^ Mq, is short enough that black holes 
can segregate from a larger volume of stars than 
in the Milky Way. In O'Leary et al. [H, we dire- 
city solved for k, and a, , using the average density 
of stars at the radius of influence and the chosen 
initial mass function of the stars. Recent obser- 
vation of the Galactic nucleus shows evidence for 
an extremely top-heavy mass function [H^. Equa- 
tion (jClOp shows that the rates are strongly dom- 
inated by the fraction of galaxies with relatively 
large values of k,. For example, galactic nuclei with 
a K, = 50% mass-fraction of compact objects (so 

that K = 20) contribute kI^^^'^ = 760 times larger 
event rate than the nominal Milky Way estimate 
with fi = 1. This explains the larger values of event 
rates in models E and F in O'Leary et al. (2009). 

• Other effects may further increase the event rates, 
which were not included in O'Leary et al. (2009). 
Keshet, Hopman, & Alexander (2009) found that 
galactic nuclei dominated by light objects can lead 
to a steeper density profile for the massive objects 
with 2 < a, < 3. Eq. (jC7p shows that the rates are 
exponentially sensitive to a,, these larger values 
lead to much higher rates for the larger mass ob- 
jects. Further, mass segregation in vector resonant 
relaxation leads to an anisotropic configuration of 
compact objects, which increases the rates by the 
square of the linear flattening of the distribution. 
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